  * Table A.5: RPM control function: s.e. clustered at the communitylevel and the household level
    set seed 101
    matrix b52 = b50[1,1..6], 0.0, 0.003, 0.26, 0.001, 0.0, 0.0, 0.0, 0.0,  0.001, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0
    mixlogit d_  dist_dw_hf_2 dist_22   nu_8d  alt2-alt4  nu_8dsq , group(choicesit) id(idid) cluster(q01_10) rand(doctor_tot  generator2009 beds p_i1 p_i2 p_i3 p_i4 ) from(b52, copy) 
    mixlpred PR0, nrep(500)
    est store mixl51
    gen esample = (e(sample)==1)
	
	
	set seed 101
    matrix b52 = b50[1,1..6], 0.0, 0.003, 0.26, 0.001, 0.0, 0.0, 0.0, 0.0,  0.001, 0.01, 0.0, 0.0, 0.0, 0.0, 0.0
    mixlogit d_  dist_dw_hf_2 dist_22   nu_8d  alt2-alt4  nu_8dsq , group(choicesit) id(idid) cluster(hh_id09) rand(doctor_tot  generator2009 beds p_i1 p_i2 p_i3 p_i4 ) from(b52, copy) 
    est store mixl51hh

